High-statistics measurement of the /3-delayed o: spectrum of Na 
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A measurement of the ^°Na /3-delayed alpha spectrum with a high-granularity set-up has allowed 
the decay scheme to be revised on several points. Three new transitions of low intensity are found 
at low a-particle energy. An i?-matrix fit of the complete spectrum gives an improved description 
of the decay and indicates feeding to the broad 2+ a-cluster state close to 9 MeV. 
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I. INTRODUCTION 

The J^ = 2+, T == 1 ground state of ^ONa de- 
cays to ^°Nc by positron emission with a half-life of 
447.9(23) ms and an available decay energy of Qec = 
13 892.5(11) keV [l|. Excited states in ^ONe populated 
in the /3 decay of ^"Na, with an excitation energy above 
4 729.84 keV, may break up to a -I- ^"^O. The a parti- 
cles originating from such transitions are referred to as 
/3-dclayed a particles and have been the subject of nu- 
merous experimental investigations. 

Through angular correlation studies, see e.g. Refs. Q, 
y] , the ratio of the Fermi (vector) and the Gamow- Teller 
(axial- vector) component has been determined for a num- 
ber of allowed transitions including the J'^ = 2+, T = 1 
isobaric analog state (IAS) at 10 273 keV. For several 
transitions to J'^ = 2+, T = states in ^°Ne, limits have 
been set on the Fermi contribution and from this the ^°Ne 
isospin mixing has also been established. 

Through measurements of ft values to individual lev- 
els in ^"Ne and comparison to the ft values in the mirror 
decay, ^°F — > ^°Nc* + 13~^ +'^e, the energy-dependence of 
the mirror asymmetry has been studied, yielding infor- 
mation on nuclear structure, meson exchange currents, 
and second-class currents, see e.g. Refs. @, Q. 

The data presented here were originally intended as 
calibration data for a precision measurement of the (3 
decay of *B @, but the quality and high statistics of the 
data reveals new features of the decay scheme of ^°Na. 
New transitions appear at low a particle energies and 
clear signatures of interference leads to a reinterpretation 
of the decay scheme at higher a-particle energies. In 
particular, we shall present evidence for the population 
of a broad resonance in -^^Ne, indicative of a a + ^^O 



cluster structure. It can be studied in a very clean way 
through (3 decay since the initial state and the transition 
operator are well understood so that the final state is 
the only unknown, see Ref. [6| for a similar study in the 
A — 12 system. 

The broad (cluster) states can play a role in the 
^^0(Q;,7)^°Ne reaction, which proceeds at a slow rate 
during astrophysical He burning due to the absence of 
natural-parity resonances inside the Gamow window. Re- 
cently Hager et al. [3] measured the total S'-factor at 
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2.26 MeV and determined the contributions from 
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direct capture transitions to the ground and first-excited 
state. Costantini et al. [3] determined S'-factor values 
for a range of energies and used an i?-matrix fit to ex- 
trapolate the S'-factor down to the energies relevant for 
astrophysics, but they left out the broad 0+ and 2"*" res- 
onances of the fourth 0+ rotational band of ^°Ne. Our 
results may allow to judge more precisely the effect on 
the S-factor of the low-energy tail of the 2^ resonance. 



II. EXPERIMENT 
A. Beam production 

The radioactive ^°Na beam was produced at the Ac- 
celerator Laboratory of University of Jyvaskyla through 
the ^^Mg(p, na)'^°Na reaction. The proton primary 
beam had an energy of 40 MeV and a typical inten- 
sity of 8-10 /zA, and the target was self-supporting 
4.3 mg/cm^ natural Mg. The activity was extracted 
with the Ion Guide Isotope Separation On-Line (IGISOL) 
technique [9|, using a light- ion fusion ion guide [lO[ and 
the mass separated A/q = 20 beam was implanted in a 
carbon foil of thickness 26 ± 2 /ig/cm^. The acceleration 
voltage was 20 kV resulting in the ions, on average, being 
implanted at a depth of 7.4 /ig/cm^. The average ^°Na 
implantation rate was 2.2 x 10^ ions per second. 

Data taking was split in two runs, the first lasting 
5 hours, the second 7 hours. In between, the setup was 
used for measuring the decays of ^^^Al (8 hours) and ^B 
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FIG. 1: Schematic and simplified illustration of the experi- 
mental setup (top view). 



tions were also made for the non-ionizing energy loss in 
the active volume of the detector which does not con- 
tribute to the observed signal. The calibration includes 
a 5 keV quadratic component due to a non-linearity in 
the electronical response (deduced from measurements 
with a precision pulse generator) and a 20 kcV quadratic 
component due to the changing pulse shape of physical 
particles at low energies. See Ref. Q for further details. 



(72 hours). At no point during this interval was the vac- 
uum broken. The ^B and ^"^Al measurements are dis- 
cussed elsewhere [a, [ii|. The two ^"Na runs give consis- 
tent results. 



B. Detection system 

The set-up consisted of four 60 uni thick, double sided 
silicon strip detectors (DSSSD) |12[ backed by 1.5 mm 
thick, unsegmented silicon detectors used to veto against 
/3 particles. The detectors were placed 5 cm from the 
carbon foil in a rectangular configuration with the carbon 
foil perpendicular to the beam, see Fig. [TJ The total 
solid-angle coverage was 30% and the angular resolution 
3 degrees. An important feature of the DSSSDs is the 
very thin dead layer of only 100 nm (over 95-96% of 
the active surface) which facilitates the detection of low- 
energy ions [l2[- The energy resolution of the DSSSDs 
was 25 keV (FWHM). 

The trigger thresholds of the data acquisition system 
were set as low as possible above the noise and were 
studied with an ^^^Am source in air at atmospheric pres- 
sure [5[ . The trigger efficiency was found to rise gently as 
a function of energy, increasing from 0% to 100% within 
an interval of roughly 100 kcV. The trigger thresholds, 
defined as the energy at which the efficiency reaches 50%, 
ranged from 160 to 240 keV, depending on the electronics 
channel. Low-energy particles below the trigger thresh- 
old could be detected if in coincidence with a higher- 
energy particle. Low-energy cutoffs in each ADC channel 
ranged from 60 to 230 keV. 



C. Energy calibration 

The two most intense /3-dclaycd a lines of ^°Na at 
2153.2(10) kcV and 4433.8(15) keV were used for the en- 
ergy calibration. Their energies were deduced from the 
excitation energies, 7421.9(12) keV and 10273.2(19) keV, 
of the corresponding states in ^°Ne and the a + ^^O 
threshold energy of 4729.84(1) keV given in Ref. [Hi. 
SRIM stopping power tables [ij] were used to correct for 
the energy loss of the a particles and the ^^0 ions in the 
carbon foil and the detector dead layer, taking into ac- 
count the variation in effective thickness with angle, but 
assuming a fixed implantation depth in the foil. Correc- 



III. DATA ANALYSIS AND RESULTS 

The singles- and coincidence a-spectra, summed over 
all four detectors, obtained in the second run are shown 
in Fig. [5] The ^^O recoil lines corresponding to the two 
most intense a lines (5 and 9) are clearly visible in the 
singles spectrum at 540 and 1 110 keV, i.e. one- fourth of 
the a-particlc energies. Only coincidence events which 
fulfill the requirement Ei6q < ■jEa + 75 keV are used 
to generate the coincidence spectrum. This cut greatly 
reduces the response tails. In particular, the response tail 
of the most intense a line (5) is greatly reduced, clearly 
exposing the new a lines below 1.5 MeV. Satellite peaks, 
caused by additional energy loss of a particles striking the 
aluminum grid covering 4-5% of the detector surface, are 
visible on the low-energy fiank of the two most intense a 
lines and are marked with asterisks. 



A. New low-energy lines 

This section is concerned with a lines 1-4; lines 1-3 
have not been observed previously. As we shall argue, 
lines 1 and 2 may be identified with a decays of the 
3-, T = state at 5 621.4(17) keV and the 1", T = 
state at 5 787.7(26) keV. The corresponding a-particle 
energies are 713 keV and 846 keV and the ^^O energies 
are 178 kcV and 212 keV. The two new a lines are only 
seen in DSSSD 3, the reason being that DSSSD 2, which 
is placed opposite of DSSSD 3, is the only detector with 
sufficiently low ADC cutoffs to detect the coincident ^^O 
ions. The "'^^O ions lose a substantial part of their en- 
ergy in the carbon foil and in the detector dead layer. At 
average implantation depth and 45° exit angle, the en- 
ergy loss is approximately 50 keV. In the active volume 
of the detector another 30 keV is lost to non- ionizing pro- 
cesses, thus reducing the detectable ^®0 energies to ap- 
proximately 100 and 130 keV, respectively, which is com- 
parable to the typical ADC cutoff in DSSSD 2. Below, 
we discuss the consequences this has for the low-energy 
a lines. 



1. Reduced detection efficiency and energy shift 

The implantation-depth distribution of 20 kcV ^°Na 
ions in carbon obtained from a TRIM simulation [ij] is 
shown by the solid line in Fig. [3] The large width of 
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FIG. 2: (Color online) /^-delayed a spectrum of '^"Na. The solid line (black) is the a + ^^O coincidence spectrum. The dashed 
line (red) is the singles spectrum. The '^^O recoil lines visible at 540 and 1 110 keV, corresponds to the two most intense a lines 
(5 and 9). Asterisks indicate satellite peaks due to the additional energy lost by a. particles striking the detector aluminum 
grid. 
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FIG. 3: Implantation-depth distributions obtained from 
Monte Carlo simulations showing the reduction in 
coincidence-detection efficiency and the apparent shift 
toward shallow implantation for the a lines at 713 keV and 
846 keV. 



the distribution causes a significant broadening of the 
1^0 energies (« 30 keV FWHM), with shallow implan- 
tations giving the largest ^^O energies because DSSSD 2 
is located upstream with respect to the foil. The re- 
coil shift discussed in Section IIII C 21 also contributes to 
the broadening (w 30 keV FWHM). Owing to the ADC 
cutoffs there exists a preference for shallow implantation 



and positive recoil shift, implying a larger-than-average 
energy loss and a negative recoil shift for the a particles. 

To determine the reduction in coincidence-detection ef- 
ficiency and the shift in the apparent a-particle energy 
due to the effects discussed above, we perform Monte 
Carlo simulations Q. The reduction in detection ef- 
ficiency and the preference for shaUow implantation is 
clearly seen in Fig. [31 which shows implantation-depth 
distributions predicted by the Monte Carlo simulations: 
If all implantations are considered, the solid distribution 
results. If, on the other hand, we only consider implan- 
tations where both the a particle and the ^^O ion are 
detected, the dotted and the dashed distributions result. 

In the first run. a lines 1 and 2 are observed at 
699(4) keV and 849(8) keV. The apparent shift in a en- 
ergy determined from the simulations are —14(5) keV 
and —4(4) keV. The systematic uncertainty on the en- 
ergy shifts was estimated by performing simulations with 
all ADC cutoffs raised and lowered by 10 keV and with 
the detector dead layer increased and decreased by 10%. 
Correcting the observed energies for the shifts, we ob- 
tain 713(6) keV and 853(9) keV. In the second run, a 
lines 1 and 2 are observed at 698(4) keV and 838(5) keV. 
The a energy shifts determined from the simulations are 
— 16(5) keV and —6(4) keV. Compared to the first run, 
the shifts in a energy are slightly larger. This is due 
to the accumulation of material on the foil between the 
two runs, resulting in a 20% increase in foil thickness. 
Correcting the observed energies for the bias, we obtain 
714(6) kcV and 844(9) keV. Taking the average of the 



TABLE I: Relative intensities of a lines 1-5 and 9. A*' is the 
number of coincidence events observed in the first (top) and 
second (bottom) run with the a particle detected in DSSSD 3 
and the ^®0 ion detected in DSSSD 2. e is the coincidence- 
detection efficiency determined from Monte Carlo simulations 
(arbitrarily fixed to 100% for line 9). 7 is the efficiency- 
corrected intensity relative to that of line 5 with the statisti- 
cal uncertainties given in parentheses. The average of the two 
runs is given in the next-to-last column and the systematic 
uncertainty on the intensity is given in the last column. 



Line 



iV E (%) 



/(%) 



Syst. uncert. 



23 24 0.043(9) 

25 17 0.023(5) 

15 65 0.010(3) 



37 


57 


0.0100(16) 


22 


89.7 


0.010(4) 


29 


87.7 


0.005(2) 


18 


93.9 


0.009(2) 


50 


93.8 


0.0082(12) 


2.139 X 10^ 


95.6 


100 



6.182 X 10^ 95.4 100 

„ 3.986x10^ 100 17.81(10) ,^^^,^. 

9 ^ 17.85(5) 

1.157 X lO'' 100 17.86(6) 



0.028(4) t'oTis 

0.0101(14) ±«:»?6 

0.006(2) ±0.0002 

0.0083(10) ±0.0002 
100 



±0.3 



two runs, we obtain 714(4) kcV and 847(5) kcV in excel- 
lent agreement with the literature values, 713.2(14) keV 
and 846(2) kcV. 



2. Relative intensities 

The determination of the intensities of the low-energy 
a lines relative to the most intense a line is presented 
in Table HI The number of coincidence events observed, 
iV, with the a particle detected in DSSSD 3 and the ^^O 
ion detected in DSSSD 2, is given in the second column. 
The coincidence-detection efficiency, e, determined from 
Monte Carlo simulations and arbitrarily fixed to 100% 
for line 9, is given in the third column. The efficiency- 
corrected intensity, /, relative to that of line 5, is given 
in the fourth column, and the average of the two runs is 
given in the fifth column. Finally, the systematic uncer- 
tainty on the intensity is given in the last column. The 
results of the two runs are generally in good agreement. 
Note the large systematic uncertainty on line 1. 

For lines 1 and 2, the systematic uncertainty is dom- 
inated by the uncertainty on the ADC cutoff's and the 
dead layer thickness. For lines 3-5 and 9, the systematic 
uncertainty is dominated by uncertainties in the model- 
ing of the ^^O ions as they strike the aluminum grid that 
covers 4-5% of the detector surface. 



TABLE IL Overview of the rare decay sequences inferred from 
the observation of the new low-energy lines. 



Line 



Decay sequence 



bfi (%) r^/r (%) 



^Ne(5.621) 



^O, 



< 0.007 



^°Ne(9.873) ^ 2°Ne(5.621) ^ ^"Ogs 0.028(14)" « 7" 
2°Ne(10.274) ^ 2"Ne(5.621) A ^^Og, 2.877(42)" 0.084(17)'' 

2 20Ne(5.788) A ^^Ogs 



'Ne(12.39) A "0(6.130) 



^O, 



0.0016(5) 
0.0010(3)" 



^Ne(10.274) A 2''Ne(6.725) A "Ogs 2.877(42)" 0.046(6)" 
^Ne(11.262) A 2''Ne(6.725) A "Og^ 0.205(26)" 6.5(8)"^ 



"Value from Ref. Il3l 

''Using r^ = 0.097(14) eV [H and T = 116(17) cV determined in 
the present work. 

"^The ratio of the phase space available for /3 decay and electron- 
capture decay (EC) to the 12 390 keV state suggests that this would 
be a 'mixed' transition with EC//9 = 0.16. 

"^Branching ratio inferred under the assumption that the decay 
sequence accounts for all the observed events. 



3. Origin of the low-energy lines 

For easier reference, our suggested decay sequences are 
summarized in Table |TT1 The detailed arguments are as 
follows. 

a. Lines 1 and 2. The energies of the two lowest- 
lying a lines are consistent with a decay of the 3~, T = 
state at 5 621 keV and the 1", T = state at 5 788 keV 
to the ground state of ^^O. The intensities are consis- 
tent with first-forbidden /3 transitions. However, it is 
possible that allowed /3 transitions to higher lying states 
followed by 7 decay to the 5 621 keV and 5 788 keV 
states also contribute. In fact, the IAS is known to 
decay to the 5 621 keV state with a partial 7 width of 
r..^, = 0.097(14) eV [H]. The a width of the IAS has been 
determined to be r„ = 116(20) eV [111. It follows that 
F-yj/Fa = 8.4(19) X 10""*. In comparison, the efficiency- 
corrected relative intensity determined from the present 
data is /i//g = 1.6(2)1q'7 x lO^'^, implying that decays 
through the IAS account for approximately 20-100% of 
the observed intensity of a line 1 ; the wide range reflects 
the large uncertainty on the efficiency correction. Fur- 
thermore the 3+ state at 9 870 keV is known Q to be 
fed in /3-decay with a branching ratio 0.028(14) and to 
7-decay to the 5 621 keV state with a branching ratio 
about 0.07 [13|, thereby giving a contribution of roughly 
the same magnitude as decays through the IAS. Note 
that the efficiency correction assumes a first-forbidden 
j3 transition to the 5 621 keV state. However, since the 
combined recoil caused by a /3 transition to the IAS fol- 
lowed by a 7 transition to the 5 621 keV state is similar 
in magnitude, the two decay modes should have similar 
efficiency corrections. We note that 7 transitions to the 
5 788 keV state have not been observed [lal . 



b. Line 3. The third a Une is located at Ea = 
1 220(30) keV, which corresponds to an excitation en- 
ergy of 6 260(40) keV in ^°Nc, assuming that the a decay 
proceeds to the ground state of ^^O. No state is known 
in ^°Ne at this energy. In fact, the observed width of Une 
3 is incompatible with such a hypothesis: Assuming a re- 
duced a width equal to the Wigner limit |17| , one obtains 
an s-wavc a width of only 7.9 keV due to the small rela- 
tive energy.^ Since no other particle-decay channels arc 
open, this represents an upper limit on the total width 
of the state. However, the observed width of line 3 is 
150(30) keV, well above the estimated experimental res- 
olution of sa 35 keV. Converting the energy scale from a- 
energy to excitation energy, this corresponds to a width 
of r = 190(40) keV, well above the Wigner limit. 

(3 particles were detected in the 1.5 mm thick silicon de- 
tectors placed behind the DSSSDs. The ratio of triple co- 
incidences {/3 + a + ^^0) to double coincidences (a-f-^^O) 
is shown in Fig. HI This '/3 ratio' displays a clear system- 
atic dependence on a energy except for line 3, which has 
a very small /3 ratio, consistent with zero at the 2a level. 
This behaviour is suggestive of electron capture (EC). 
We have identified one possible decay sequence: The 1+, 
13 308 keV state in ^ONe is only fed in EC decay. It a 
decays to the 3-/2+/1- states in ^^O (cf. Fig. ^ with 
partial widths of 700/19/190 in units of eV, emitting a 
particles with energies of 1960/1330/1170 in units of 
keV. The main a line at 1 960 keV would be buried under 
the tails of higher-lying a lines. However, the observed 
intensity of line 3 implies a log ft value of 2.3, an ex- 
tremely unlikely possibility. We have searched for other 
potential candidates among the known states in ^"Ne, 
but have failed to find any. The origin of line 3 thus 
remains mysterious. 

We note, however, that the sd-shell model calculations 
of Ref. [Hi predict both a 3+ state at 11400 keV and 
a 1+ state at 12 200 keV to be populated in /3 decay. 
None of these states have been observed in /3 decay yet, 
but (3 decay to the 1+ state with the calculated strength 
followed by a decay to the 3~ state in -"^^O would explain 
the properties of line 3, provided the l"*" state lies at 
the slightly higher energy of 12 390 keV. At this elevated 
energy (3 decay will still be favored over EC decay with a 
phase-space ratio of EC//3 = 0.16. It would thus appear 
that the proposed decay sequence is inconsistent with the 
anomalous low /3 ratio of line 3. However, the parameter 
driving the systematic behavior seen in Fig. 4 is the (3 
energy, not the a energy. In the proposed decay sequence 
the /3 energy is much smaller than the a energy would 
indicate because the decay feeds an excited state in ^^O; 
if the decay were to feed the ground state the a energy 
would be 6 130 kcV. Indeed, the /3 ratio of line 3 is in good 
agreement with the observed systematics at 6 130 keV. 
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FIG. 4: Ratio of triple coincidences [13 + a + ^®0) to double 
coincidences (0:+^ O). The dashed curve is only meant to 
guide the eye. The error bars show the statistical uncertainty. 



We note that the a threshold in ^^O is only 1 MeV 
above the 3~ state and that, in principle, population of 
unbound states in ^^O is possible in the decay of ^'^Na. 
We have searched for such /3-delayed two-a events, but 
found no significant signal above background. 

c. Line 4- In Ref. [3| it was suggested that line 4 
should be attributed to the a decay of the 0+ , T = state 
in ^°Ne at 6 725(5) keV. The excitation energy deduced 
from the present data is 6 716(6) keV, thus supporting 
this interpretation. The intensity is a factor of two be- 
low the intensity obtained in Ref. Q. We arc unable to 
account for this discrepancy. As argued in Ref. Q , the di- 
rect population of a 0^ state requires a second-forbidden 
transition, which results in a much smaller intensity than 
the intensity actually measured. Two alternative expla- 
nations were offered: (i) a /3 branch feeding a hypothet- 
ical state at 12 758 keV, which a decays to the 0+ first- 
excited state in ^^O at 6 049 keV, and (ii) a weak /3 branch 
feeding a hypothetical state around 11 000 keV in ^°Ne, 
which decays mainly by 7-ray emission to the 6 725 keV 
state. It was already noted in Ref. [3| that alternative 
(i) is unlikely to be the correct explanation because the 
f3 branch would have a very small log/i value of 3.28. 
We can fully rule out alternative (i) : Since the 6 049 kcV 
state in ^^O decays by e'''e~ emission, the (3 ratio should 
be significantly higher (roughly a factor of 3) for a line 
4 compared to all other a lines. Fig. 0] clearly shows 
that this is not the case. We cannot confirm nor rule 
out alternative (ii). Instead, we offer a third and fourth 
alternative: (iii) a weak 7 branch from the IAS to the 
6 725 kcV state and (iv) a weak 7 branch from the 1+, 
T =1 state at 11 262.3(19) keV to the 6 725 keV state. 

We first consider alternative (iii): Based on the ob- 
served intensities, we can easily estimate the width of 



the proposed 7 branch, 



Tn 



h 



0.054(11) cV 



(1) 



where I^ and /g are the efficiency-corrected relative in- 
tensities given in Table U and Fq. — 116(20) eV is the a 
width of the IAS [1^1 ■ The least intense 7 branch ob- 
served by Fifield et al. was that to the 4 967 keV with a 
partial width of 0.060(8) eV [H. It is difficult to judge 
from their Fig. 2 whether a 0.054(11) eV branch to the 
6 725 keV state could have gone unnoticed, but it ap- 
pears possible, especially taking into account that the 
6 725 keV state is fairly wide (F = 19.0(9) keV). 

Next, we consider alternative (iv): The 11262 keV 
state is fed with a /3-branching ratio of 6/3(11.262) = 
0.205(26)% corresponding to a log/t value of 3.72 Q. 
a decay from the 11 262 keV state is parity forbidden^ so 
it almost exclusively decays by 7 emission. It is known 
to decay to the ground state and the first-excited state 
with relative intensities of 84(5)% and 16(5)%, respec- 
tively [201 , and the partial width to the ground state has 
been determined to be 11(2) eV [21[. Based on the ob- 
served intensity of line 4, we can estimate the width of 
the proposed 7 branch to be 



bp{JA2l) 
^^" "^6^11.262) 



= 0.086 eV 



(2) 



where I4, is the efficiency-corrected relative intensity 
given in Table H F^ = 13.3(2) eV is the total width 
of the 11262 keV state and 6/3(7.421) = 15.96(22)% is 
the /3-branching ratio to the 7421 keV state. The de- 
duced width corresponds to an Ml transition strength of 
0.043 W.U., which is a reasonable value. Furthermore, 
the branch is far smaller than the two known branches, 
so it could easily have been missed in the study of Berg 
et al. [20|; unfortunately, their 7 spectrum only extends 
down to 6 MeV. 



B. Energies and relative intensities 

The literature data on the energies and relative in- 
tensities of the /3-delayed a lines of ^"Na are from 
Refs. d,!!,!!!-!!!. in Tablejml we compare the present 
data to that of Refs. 0, [23, [23l- The lines are numbered 
as in our spectrum and letters are used to denote lines 
quoted in other papers, but not seen by us (as argued 
below most of these are due to misinterpretations). Line 
D, seen by Ref. [23| but not by Ref. jj], was not seen 
in the present study despite the higher statistics. The 
present data on the energies of lines 4-13 is consistent 
with the previous data. The energies of the two most 



intense a lines (5 and 9) used for the energy calibration, 
differ from the previous energies by a few keV. This is 
due to the updated value for the a 4- ^^O threshold in 
20Ne Up. 

The intensity of line 9 agrees with the intensity re- 
ported in Ref. (2J^^but is slightly larger than the intensity 
reported in Ref. Q. An older study [J] gives 17.39(11)% 
in agreeement with Ref. 13|. 

Only the low-energy lines and the two most intense 
lines are included for the present data. The less intense 
lines at higher energy can be influenced in position as well 
as intensity by interference effects so we shall give our 
results below, when we discuss the i?-matrix fit (Table 
lYl). 



C. _R-Matrix Analysis 

We turn now to a detailed analysis of the /3-delayed a- 
spectrum above 2 MeV. This part of the spectrum con- 
tains several resonances as well as indications of inter- 
ference effects at lower intensity levels. The _R-matrix 
formalism provides a suitable framework to incorporate 
such effects in a fitting routine, but reliable results can 
only be extracted once the experimental response func- 
tion is fully understood. 



1. R-matrix parametrization 

The i?-matrix parametrization allows extraction of 
energy-, width- and /3-strength parameters for the res- 
onance levels in question. The parametrization offers at 
the same time a more proper interpretation of the excita- 
tion energy distribution compared to earlier studies. We 
use the same procedure as in Ref. |5| . The probability for 
/3 decay to the excitation energy E, followed by breakup 
into the channel c, is given by 

w, (E) = C^fp [E] P, {E - Sthrcs) 



x=F,GT Xfi 



(3) 



where the alternative parametrization of the level matrix 
Axfj^ is used, as given in Eq. 33 of Ref. |26[. In Eq. ([3]) the 
boundary parameter Be is automatically set equal to the 
shift function Sc{E\) for all levels. This ensures that the 
fitted parameters for the energy. Ex , the /3 strength, g\x , 
and reduced a width, 7^^ , becomes the so-called observed 
ones [20]. The notation x =F, GT refers to Fermi and 
Gamow- Teller decays respectively, fp is the integrated 
phasespace available to the leptons, C is a normalization 
constant and Pc is the penetration function with -Ethrcs 
being the a-threshold energy. The number of counts as 
a function of E is then 



^ The parity-forbidden a-docay width has been determined to be 
42(20) X 10-'5 eV [11. 



NiE) 



In 2 



[E). 



(4) 



TABLE III: Comparison of present energies and relative intensities to previous data. 



Line 




E (keV) 






h (%) 








Present 


Ref. [31 


Ref. [23| 


Present 


Ref. \S} 


Ref. [23] 


Ref. [24] 


1 


714(4) 




~ 780" 


0.028(4)(ttS) 




~ 1.4" 




2 


847(5) 






0.0101(14)(lf6) 








3 


1220(30) 






0.006(2) 








4 


1589(5) 


1 580(40) 




0.0083(10) (2) 


0.020(4) 






5 


2153.2(10)'' 


2 150.4' 


2 150.4' 


100 


100 


100 


100 


6 




2 479.6(21) 


2 483.5(25) 




3.7(4) 


3.91(24) 


3.7(3) 


A 




2 659(7) 


2 756(5) 




0.074(6) 


0.072 




7 




3 570(25) 


3 325(13) 




0.39(4) 


0.052 




8 




3 799(3) 


3 803.0(25) 




1.510(27) 


1.55(5) 


1.6(2) 


9 


4433.8(15)' 


4 432.2' 


4 432.2' 


17.85(5)(30) 


17.31(9) 


17.83(27) 


16.3(3) 


10 




4 675(3) 


4 674.6(21) 




0.553(15) 


0.478(24) 


0.45(6) 


11 




4 885(3) 


4 884.4(25) 




1.09(3) 


1.15(4) 


0.91(8) 


B 




4 966(7) 


4 930(6) 




0.075(9) 


0.128 




C 




5 106(7) 


5 222(47) 




0.055(7) 


0.044 




12 




5 249(4) 


5 253.5(23) 




0.165(11) 


0.130(9) 




13 




5 698(6) 


5 691(4) 




0.010(2) 


0.012(2) 




D 






5 896(6) 






0.002 





"The evidence presented in Ref. 
group is meager. 
'Used for energy calibration. 



I23II for the observation of this a 



where N is the total number of counts and Ti /2 is the 
half life of 20Na [H. 

The observed a widths are extracted as 



n 



Ec2Pc7L 



i + Ec7LMb=i.. 



(5) 



and the procedure for determining the /3-decay matrix 
elements is similar to the one in Ref. [23] 



B 






E, 



thrc 






dE. 



(6) 



Here B = 6144.2(1.6) s ^ and Q is the ^-decay Q- 

2 
value. For B\gt a factor of ^- enters in the denomina- 

tor. When calculating absolute Bp/^T values the known 
decay branch (BRq, — 20.05 ±0.36 Q) to a unbound lev- 
els in ^°Ne is used. Integration of eq. ([3]) (without C^) 
over the full energy spectrum then yields the normalisa- 
tion constant 



C 



BR„ In 2 



[T,/,J^w{E)dE^ 



1/2 



(7) 



In our analysis only the exit channel populating ^^O in 
the ground state has been included. Levels 11, 12 and 13 



are above the threshold for decay to exited states in ^^O. 
Their contribution to the reduced a- widths have been in- 
vestigated by performing a two-channel fit including the 
first exited 0+ state in ^^O at 6 049.4(1) keV. From this 
test their contribution was found to be no higher than 
10"'^ VkeV which justifies the single-channel analysis. 



2. Recoil Broadening 



If the 20Ne* 



16 



O breakup occurs at rest, the 



energies of the a particle and the ^^O ion are fixed by 
energy and momentum conservation. The available en- 
ergy is shared as EieQ/E^ = ma/mi6Q = 0.25. If the 
breakup is preceeded by a /3 decay, however, the energies 
of the a particle and the ^^O ion are smeared out due to 
the recoil motion of the ^°Nc* nucleus. The recoil motion 
also causes a small systematic shift of the mean energies 
which, however, may be safely neglected here. For a pure 
Fermi transition, the resulting energy distribution can be 
approximated by [29|, 



p{x) = «^"^ 



-(l-x4) , -1 < X < 1 
, |a:| > 1 



(8) 



while, for a pure Gamow- Teller transition and the spin 
sequence 2+ -)► 2+ -> 0+, 



p{x) 



16T„, 



-(l-2a;2+a;4) -1 < x < 1 
\x\ > 1 



(9) 
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FIG. 5: Recoil broadening, (a) Maximum recoil shift as a 
function of the excitation energy in Ne. (b) Fermi (F) and 
Gamow- Teller (GT) recoil broadening distributions. 



with X — 5E/Tjnax where SE = Ea — (Ea) is the shift 
relative to the mean a-particle energy and T^ax is the 
maximum shift, given by 



T„ 



M 



2Qm^c^{W^ - 1 



il/c2 



Q 



Mc^ - Q 



1/2 



where 

20 



me 



(10) 
and AI are the electron, a-particle, 
and 20Ne masses; Q = E^ - A 729.84 keV and Wq = 
{Eq — Ex)/meC^ where E^ is the excitation energy in 
^°Ne and Eq = 13 376 keV is the maximum total /3 en- 
ergy for decays to the ground state of ^°Na. In Fig. [S] (a) 
the maximum recoil shift, T,„ax, is shown as a function of 
the excitation energy in ^"Ne, E^- The inset (b) shows 
the different shapes of the recoil broadening distribution 
for Fermi (F) and Gamow- Teller (GT) transitions. 



Response Function 



is therefore well understood; this allows for a detailed 
study of the modification of the shape due to experimen- 
tal effects. 

We adopt a parameterization similar to the one used in 
Ref. [201 to describe the experimental response function. 
It consists of a Gaussian and a Gaussian folded through 
a low-energy exponential tail: 



^{E^,E) = 



Ai {E- EoY 

exp ' 




2a2 



X erfc 



E-Eo + o-yA 
V2a 



(11) 

Here Eq and E are the nominal and observed energies, A 
the exponential decay "length" and erfc the complement 
of the incomplete error function. The normalization con- 
stants are Ai = 1/(1 + r) and A2 = r/(l + r) where r 
gives the relative size of the second term compared with 
the first one. We find A = 25.5 keV and r = 0.085. In 
addition we allow for contributions from low energy tails 
of higher lying Gamow- Teller fed resonances. 

Ions striking the aluminum grid covering g = 4.9% 
of the detector surface experience additional energy loss 
compared to other ions, thus giving rise to a satellite 
peak at slightly lower energy. (The fit parameters are 
fairly strongly correlated so that almost equivalent fits 
can be obtained with different sets of parameters; we 
quote here a consistent set corresponding to the fit op- 
timum.) The satellite is wider than the main peak due 
to the larger variation in effective thickness with angle. 
The final complete response function therefore reads: 



^iEo,E) ={l-g)xJj{a;Eo,E) 

+ 9i^{(^g;Eo- Eg,E) , 



(12) 



where Eg is the mean energy loss in the grid, a = 
10.2 keV the Gaussian width of the main peak and 
as ~ 25.2 keV the Gaussian width of the satellite peak. 
The shape of the IAS a line, including experimental 
effects, is found by convolution of the response function, 
Eq. p^ . with the recoil broadening distribution, Eq. ([S]), 



The response function describes the distribution of en- 
ergies measured from a perfectly monochromatic source 
due to experimental effects. We describe now how the re- 
sponse function is extracted from the experimental data. 

The second most intense line (9) in the /3-delayed a 
spectrum of ^"Na is found at 4 434 keV and results from 
transitions to the IAS in ^"Nc at 10 273.2 keV, the width 
of which is less than 0.3 keV [I3. The transition to the 
IAS is known Q to be a mixed Fermi and Gamow- Teller 
type transition. However, the Gamow- Teller component 
is rather small so the broadening effect is well described 
by Eq. dH) . Owing to the special structure and the T = 1 
nature of the IAS, only little interference is expected with 
neighboring 2+ states which could potentially distort the 
shape of the a line. The physical shape of the IAS a line 



dN 
dE 



j ^{E^,E)p{x{E^))AE^ 



(13) 



Fig. [6] shows the line shape measured in DSSSD 1 with 
the recoil broadening distribution, the response function, 
and the best fit using Eq. (|13|) superimposed. 

The response function, "^{Eq^E), is hereby obtained 
at 4434 keV. In order to determine the response func- 
tion for a general a-particle energy, we (1) assume that 
the exponential tails are independent of energy; (2) con- 
vert the energy loss in the aluminum grid, E'g, into an 
equivalent thickness (0.60 /xm) which we use to calculate 
the energy loss in the grid for other a-particle energies; 
and finally (3) assume the Gaussian width of the satellite 
peak, Ug, to be proportional to Eg. 



> 






^10 

o 

o 


E 


Data 


in 


- 


Fit 


o 


- 


— Gamow-Teller tail 
Recoil 




10-^ - 



10 



4.1 4.2 4.3 



4.4 4.5 4.6 

a energy (MeV) 



FIG. 6: (Color online) Singles spectrum between 4 and 
4.6 MeV measured in DSSSD 1. The short-dash-dotted 
(green) curve shows the recoil broadening distribution, 
Eq. (Is}. The long-dash-dotted (blue) line shows the contri- 
bution from GT transitions to the low-energy tails of higher- 
lying states, modeled as an exponential function. The dashed 
(red) curve shows the best fit to the data using Eq. (|13|l . 



4. Synthesis and Fitting 

The ft values obtained in Ref. Q suggest that the 
transitions to the observed resonanees are first forbidden 
or allowed, and hence the spin- parity could be 1~, 2+ or 
3~. (The observation of a decay from these states ex- 
cludes unnatural spin-parity.) We shall attempt first the 
simplest possible fit based exclusively on allowed transi- 
tions to resonances seen earlier in reaction experiments. 
All resonances will then have spin-parity 2+ and inter- 
ference occurs naturally. 

The delayed particle spectrum recorded in DSSSDl 
during the second ^"Na run has been fitted using the 
maximum likelihood estimator |3l| for the excitation- 
energy region 7.00-12.00 MeV. Figure [7] displays both 
the fitted i?-matrix function (eq. 2]) and the final fit 
curve obtained after convolution with the response func- 
tion and the recoil broadening distribution. The devia- 
tion between the data and the fit is largest at the satel- 
lite shoulder of the most intense peak and dominates the 
rather high x^/d.o.f. = 2.21. This deviation is caused 
by systematic errors in the propagation of the response 
function to other energies. The inability of the fit to de- 
scribe the few events around 9.8 MeV is due to inaccurate 
modeling of the low-energy response tail of the IAS. The 
fit results are listed and compared to literature values in 
Table IVl 

As starting point in the fitting procedure the fairly 
well known level energies were fixed whereas the more 



uncertain reduced a-widths and some /3-strength param- 
eters were left free. Eventually all parameters were var- 
ied except for the IAS reduced a- width, 79, which in the 
literature has the upper limit Ftot < 0.3 keV. Due to 
lack of sensitivity to such a small width 79 was kept at 
-6.09 X 10^3 VMcV corresponding to Ttot = 0.200 keV 
throughout the minimization. To reproduce the data it 
was necessary to include a broad level at 8 767(23) keV; 
we shall comment on this below. 

The uncertainties on the i?- matrix parameters a (Ex), 
(7 (jx) and cr (5a,f/gt) were calculated by adding their 
individual errors in quadrature. The error originating 
from the response function was estimated by comparing 
fits using two different response functions: One that in- 
cludes the low energy tail of the IAS and one that does 
not. The detailed error budget is as follows. 

Energy parameters (Ex): Errors on Ex from the 
energy calibration range from 1 .4 keV to 3.0 keV while the 
statistical fit errors range from 2 x 10""^ keV to 4.5 keV. 
The response function errors on Ex are in general small 
but relatively large for levels 7 and 13. 

Reduced width parameters (7a): The response er- 
rors dominate since the reduced widths are highly sen- 
sitive to the detector response. The errors are around 
2-4 X 10~'^ Vke V fo r all levels except level 7 for which 
it is 2 X 10^^ VkeV. The statistical errors are at the 
10^'^ VkeV level. The uncertainty on Ftot, a is found di- 
rectly by combining cr (7a) and eq. ([S]). 

/3-strength parameters (gA,G/FT): The statistical 
errors dominate giving contributions around 10"'^ for all 
levels apart from levels 6, 12 and 13 where they are one 
order of magnitude larger. The uncertainty on Bp/cT is 
obtained through propagation of errors in eq. ([5]) where 
a (BRa) is the dominating source of uncertainty. 



5. Discussion 

Our observed level parameters in Table ITVl will now be 
discussed and compared to the most recent ^"Ne evalu- 
ation [13| . Compared to earlier work we propose a new 
interpretation of the /3-delayed a-spectrum from ^°Na. 

The energy positions of levels 5, 9, 10, 11, 12 and 13 
found in this work are in agreement with the literature 
values within one standard deviation. Levels 6, 7 and 
8 differ slightly more, we note that their energy region 
contains clear interference effects. 

The literature widths in Table I IVl are total widths, 
many are quite uncertain so that a quantitative compar- 
ison to our values is not possible. For levels 5, 7, 8, 
10, 11 and 12 our widths are more precise than the lit- 
erature ones. The width of level 5 differs significantly 
from the one i n |l3| but agrees well with a newer value 
of 9.0(13) keV jg. In the case of level 6 the value of 
180(80) eV is too small to be trusted and our poor de- 
termination for level 13 is attributed to the low level of 
statistics in this part of the experimental spectrum. 

The calculated -Bp/GT values are compared to those 
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FIG. 7: Best fit to the experimental data. The graph shows the fitted R-matrix expression (eq. |4]) (dashed), the convoluted 
expression (long-dashed) and the data (solid). 

TABLE IV: Results expressed as observed _R-matrix parameters E\, Fq and Bp/Qx with comparison to literature values. 







This work 




Compilation [13] 


Reference [3| 


Level 


E[MeV] 


Fc[keVl 


Bf/gt 


E[MeV] 


Ftotai[keV] 


Bf/gt 


5 


7.4227(15) 


10.0(5) 


0.246(4) 


7.4219(12) 


15.1(7) 


0.240(8) 


6 


7.8301(14) 


0.18(8)" 


0.017(7) 


7.8334(15) 


2 


0.0126(5) 


7 


8.767(23) 


686(56) 


0.00230(7) 


9.00(18) 


800 


0.008(12)' 


8 


9.4927(19) 


35(3) 


0.0329(8) 


9.483(3) 


29(15) 


0.032(2) 


9 


10.274(2) 


0.2 


1.95(3) 


10.2732(19) 


< 0.300 


2.06(6) 


10 


10.587(2) 


34(2) 


0.069(2) 


10.584(5) 


24 


0.067(5) 


11 


10.842(2) 


16.4(5) 


0.258(5) 


10.843(4) 


13 


0.224(15) 


12 


11.331(3) 


85(6) 


0.127(4) 


11.320(9) 


40(10) 


0.101(14) 


13 


11.89(5) 


252(70) 


0.11(10) 


11.885(7) 


46 


0.039(19) 



"The systematic uncertainty is larger than the quoted fit uncer- 
tainty. 
''In [il the level is located at 9.196(30) MeV. 



from Clifford et a/., which were obtained by conversion 
from the reported log/t values. Agreement within one 
standard deviation is only obtained for levels 5,6,8 and 
10 albeit level 9 nearly agrees. The large difference for 
level 13 is again ascribed to the few events collected in 
that region while the disagreement for levels 7, 9, 11 and 
12 largely has to do with our reinterpretation of the a- 
spectrum. 

The analysis performed by Clifford et al. does not 
take interference effects into account. The events be- 
tween level 11 and 12 in their interpretation therefore 
suggest the presence of two new levels in ^"Ne (level 10 



and 11 in figure 4 in [3|), levels that have not been seen in 
other reactions. Our analysis, which treats interference 
effects properly, shows that constructive interference be- 
tween levels 11 and 12 naturally reproduces the spectral 
shape in this region. 

Clifford et al. explained the characteristic structure 
around 8-9 MeV by postulating the existence of two 
new levels in ^ONe at 8 058keV and 9 196 keV, but they 
do mention that interference between the 9 196 keV level 
and other levels could explain the 8-9 MeV structure. 
Our analysis does not include the 8 058keV level (only 
reported by Clifford et al.). The interference effects are 
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FIG. 8: The _R-matrix expression (eq. |4| including only levels 
5 and 7 plottet against data. This is not a fit, but the best-fit 
parameters for the two levels are used. Shown is the R-matrix 
expression (eq. [4]) (dashed), the convoluted expression (long- 
dashed) and the data (solid). 



here very important and are mainly due to destructive 
interference between level 5 and 7. This is shown explic- 
itly in figure [8] where eq. (HJ including only levels 5 and 
7 is plotted. 

Clifford et al. briefly mention the possibility that the 
9 196 keV level identified by them could be identical to a 
known broad level at 8 800 keV. Indeed a level with 2+ 
had been observed at « 8 800 keV both in ^^0{a, a) [H] 
and ^^F(^He,d) [3J] reaction experiments. Later Burlein 
et al. [3J] measured the ^°Ne(7r, tt') reactions, reported 
E = 9 000(200) keV and T = 800 keV for this level and 
identified it as the 2"*" member of the 0^ rotational band 
in 20Ne. The present resuh, E = 8 770(20) keV and T = 
686(56) keV, is the most precise observation of the 2^ 
level and the first conclusive observation in /3-decay. 

The pronounced ^^O -|- a cluster structure of ^°Ne has 
been verified both experimentally [35-37] and theoreti- 
cally [3^ and several rotational bands containing reso- 
nances with very large a widths have been identified. 
The properties determined for the 8 770 keV resonance 
in our work clearly indicate cluster structure. Its large 
a-width points to a substantial overlap with the exit 
channel wave function indicating a large ^^O -f- a com- 
ponent. More quantitatively 7^ = 0.320 MeV is close to 
the Wigner limit 7.^ of 0.344 MeV, evaluated with a ra- 
dius of ro = 1.5 (16^/^ -I- 4^/^) fm. The cluster structure 
is consistent with the small Bqt value determined in this 
work; an allowed Gamow- Teller transition to a pure ^^O 
+ a cluster state would be Pauli blocked. 

As a final remark we note that the position for level 
7 agrees fairly well with the position found in a theo- 
retical study J38| for the 2^ state, thus supporting our 



assignment. 



6. IAS properties 

Our results give, when combined with earlier experi- 
ments, improved information on the IAS. First we con- 
sider the fraction of delayed a-particlcs coming from the 
IAS. This is best estimated from the singles spectrum, 
Fig. [21 by including the intensity above 1.5 MeV. This 
lower limit is varied by 100 keV to estimate the system- 
atic uncertainty. The DSSSD detectors have a singles 
a particle efficiency that varies slightly with energy and 
has been determined in Ref. [5[ , the data are corrected for 
this effect. Each DSSSD gives an independent value for 
the ratio. The four values are consistent and lead to a 
final value of IfAsl^tot = 0.1391(5). Compared to the 
previous value of 0.1386(33) [3| the precision is improved 
significantly. 

Our second result is the total j3-a strength for the 
IAS of 1.95(3) given in Table |TVl The /3-a strength is 
a product of the /3 strength and the branching to the 
a decay channel and must be corrected with the fac- 
tor F^^'^/rf^'^. This factor can be estimated in the 
following way. The total Fermi strength is 2(1 — 5c) 
where the correction factor according to theory 0] is 
0.014(5). Combining this with the ratio of Gamow- 
TcUer to Fermi strength, also measured in Q, gives a 
value for the total beta strength to the IAS that leads to 

ri-"^/rL1^ = 0.92(3). 

Alternatively, one may determine this ratio by mea- 
suring the 7-decays of the IAS. Several branches have 
been seen [l^ , the main one goes to the 1 630 keV level 
and constitute 88.9(5)% of all IAS 7-decays. (This num- 
ber includes only observed transitions; as yet unobserved 
branches, such as the one to the 6 725 keV level discussed 
above, are not expected to change the value significantly.) 
Two relative intensities, that of the IAS transition to the 
1 630 keV level relative to the total number of 1 630 keV 
7-rays [33 and that of the total number of a-particles 
relative to the 1 630 keV 7-ray [j] , now suffice to extract 
^'o^^ l^iot^ = 0.961(5), a value more precise than the 
above one. 

Finally, one may include the measurements of a- 
radiative capture through the IAS to the 1 630 keV 
state [1^, [ia, [io], which give an average value of 



5r^^^i ggF^^^/F^^^^ of 19.5(1.5) eV [Ij, to put the rel- 
ative measurements on an absolute scale. This gives a 
Fj^'^ of 116(17) eV, which in hindsi ght justifies our use 
above of the 116(20) eV deduced in [la|. The difference 
to the value of 200 eV used in the fit (Table |IV| will not 
affect the fit quality. 



IV. SUMMARY AND CONCLUSION 

The serendipitous improvements to the /3a part of the 
decay scheme of ^"Na (see Fig. ^ presented here fall in 
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two groups. 

By measuring the a particle and the ^^O recoil ion in 
coincidence, we were able to identify three new /3-delayed 
a groups below 1.5 MeV that, in the singles spectrum, 
are buried under the up to five orders of magnitude more 
intense ^^O recoil groups. The two lowest groups could be 
due to first forbidden transition (they arise from known 
negative parity states in ^°Nc), but at least the lowest one 
has an important contribution from (3^a decays through 
higher lying ^°Ne states. The origin of the third group 
is less certain, but it could be due to a combined EC- 
alpha-gamma and bcta-alpha-gamma decay through ex- 
cited states in ^^O. A fomth low-energy group is also due 
to j3^a decays. 

The high level of statistics gathered in the present ex- 
periment allowed for an i?-matrix analysis of the inter- 
ference features seen in the /3-dclaycd a spectrum around 
3 MeV and 5 MeV. The spectrum turns out to be describ- 
able in terms of levels that all have been seen in reaction 
experiments. Compared to earlier /3-decay studies we re- 
move three suggested levels, but confirm conclusively /3 
feeding to the broad state close to 9 MeV excitation en- 
ergy in ^°Ne. The large width of this state suggests a 
pronounced a + ^^0 cluster structure; our analysis gives 
improved values for its energy and width. 



Most remaining questions on the ^"Na decay are minor 
and involve 7-decay. A dedicated study where a-7 coinci- 
dences were recorded may clarify the situation concerning 
the observed new low-energy a groups. Also of interest 
could be a search for population of unbound states in ^^O 
by detections of two-alpha events with a setup similar to 
the one used in the present work. 
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FIG. 9: Decay scheme for ■^°Na. Details of several rare decays are given in Table |TT1 Gamma rays in ^"Ne are not marked. 



